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An effective field tfieory is used to give a model-independent description of Comp- 
ton scattering at energies comparable to the pion mass. The amplitudes for scat- 
tering on the proton and the deuteron, calculated to fourth order in small mo- 
menta in chiral perturbation theory, contain four undetermined parameters that 
are in one-to-one correspondence with the nucleon polarizabilities. These polariz- 
abilities are extracted from fits to data on elastic photon scattering on hydrogen 
and deuterium. For the proton we find: a.p = (12.1 ± 1.1 (stat.))l[JJ (theory) and 
j3p = (3.4 lb 1.1 (stat.))^g} (theory), both in units of 10~^ fm'^. For the isoscalar 
polarizabilities we obtain: = (13.0 ± 1.9 (stat.)) jif'g (theory) (in the same units) 
while 13m is consistent with zero within sizable error bars. 



Electromagnetic polarizabilities are a fundamental property of any composite object. For 
example, atomic polarizabilities contain information about the charge and current distri- 
butions that result from the interactions of the protons, neutrons, and electrons inside the 
atom. Protons and neutrons are, in turn, complex objects composed of quarks and glu- 
ons, with interactions governed by QCD. It has long been hoped that neutron and proton 
polarizabilities will give important information about the strong-interaction dynamics of 
QCD. For example, in a simple quark-model picture these polarizabilities contain averaged 
information about the charge and current distribution produced by the quarks inside the 



nucleons In this paper we use an effective field theory (EFT) of QCD to extract both 
proton and neutron polarizabilities in a consistent and systematic manner from Compton 
scattering data — the first EFT extraction of all these important quantities within the same 
framework ^. The background to this work will be found primarily in two papers, Ref. 
for the proton, and Ref. Jj] for the deuteron. Further details of interest to the specialist will 
be published elsewhere p. 

In an atomic or molecular system the polarizabilities are measured with static fields. 
Nuclear polarizabilities can analogously be determined by the scattering of long-wavelength 
photons. Experimental facilities which accurately measure the energy of a photon beam 
using photon tagging have made possible a new generation of experiments which probe the 
low-energy structure of nucleons and nuclei. In particular, photon tagging can be used to 
measure Compton scattering on weakly-bound systems, since it facilitates the separation of 
elastic and inelastic cross sections. At sufficiently low incoming (outgoing) photon energy a; 
{uj') and momentum k {k'), the spin-averaged Compton scattering amplitude for any nucleus 
is, in the nuclear rest frame: 

/ z'^e'^ \ - - 

T = e'-e\ +ATTaujuj'] + Att/3 e' x k' ■ ex k + . . . , (1) 
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where e and e' are the polarization vectors of the initial and final-state photons. The first 
term in this series is a consequence of gauge invariance, and is the Thomson limit for low- 
energy scattering on a target of mass and charge Ze. The coefficients of the second 
and third terms are the target electric and magnetic polarizabilities, a and (3, respectively. 
The polarizabilities can be separated by the angular dependence: for example, at forward 
(backward) angles the amplitude depends only on a + (3 {a — (3). Other terms, represented 
by ". . .", include higher powers of energy and momentum and relativistic corrections 0]. 

Hydrogen targets are used to determine proton polarizabilities ap and /9p 0i- By contrast, 
the absence of dense, stable, free neutron targets requires that the neutron polarizabilities 
a„ and /3„ be extracted from scattering on deuterium (or other nuclear) targets. Data exist 
for coherent — >■ 'jd from 49 to 95 MeV 0, 10, 11 1, and for quasi- free 'yd 'jpn from 200 



to 400 MeV [12] . The coherent process is sensitive to the isoscalar nucleon polarizabilities — 
aN = {ctp + an)/2, (3n = {Pp + Pn)/'^ — via interference with the larger Thomson term. The 
extraction of these polarizabilities from data requires a consistent theoretical framework 
that clearly separates nucleon properties from nuclear effects. In the long-wavelength limit 



pertinent to polarizabilities EFT provides a model-independent way to do exactly this 

Q,li3- 

The EFT of QCD relevant to the low-energy interactions of a single nucleon with any 
number of pions and photons is known as chiral perturbation theory (^PT). The Lagrangian 
is constrained only by approximate chiral symmetry and the relevant space-time symmetries. 
S'-matrix elements can be expressed as a simultaneous expansion in powers of momenta and 
the pion mass (collectively denoted by Q) over the characteristic scale of physics not included 



explicitly in the EFT. Many processes have 
and it has proven remarkably successful [l4|. 



_j3een computed in this EFT to nontrivial orders 
In this, as in any EFT, detailed information 
about short-distance physics is absent. The short-distance physics relevant to low-energy 
processes appears in the theory as constants whose determination lies outside the scope of 
the EFT itself. In the purest form of the theory they are determined by fitting experimental 
observables. In many cases the sparsity of low-energy single-nucleon data makes such a 
determination problematic — and for free neutrons data is non-existent. 

When the amplitude for unpolarized Compton scattering is expanded in powers of Q, 
we obtain Eq. with a and (3 given as functions of the EFT parameters. To 0{Q^), 
no parameters appear apart from those fit in other processes, so predictions can be made. 
At this order the proton and neutron polarizabilities are given by pion-loop effects Q|: 
ttp = a„ = iO/5p = 10/3„ = 5e2^y3847rV|m^ = 12.2 x IQ-^fml Here ^ 1.26 is the 
axial coupling of the nucleon and /tt — 93 MeV is the pion decay constant. At 0{Q^) there 
are new long-range contributions to these polarizabilities. Four new parameters also appear 
which encode contributions of short- distance physics to the spin-independent polarizabilities. 
Thus, minimally, one needs four pieces of experimental data to fix these four short-distance 
contributions, but once they are fixed xPT makes model-independent predictions for Comp- 
ton scattering on protons and neutrons. 

The amplitude for Compton scattering on the nucleon has been computed to 0{Q^) in 
Ref. P]. The shifts of ap and Pp from their 0{Q^) values were not fitted directly to the 
data in that work. Instead the Particle Data Group values for the polarizabilities were 
used |17]. These were originally extracted using a dispersion-theoretic approach that incor- 
porates model-motivated assumptions regarding the asymptotic behavior of the amplitude. 
The differential cross sections which result from this procedure are in good agreement with 
the low-energy data 
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FIG. 1: Results of the OiQ'^) EFT best fit (solid line) to the differential cross sections for Compton 

WW 

scattering on the proton at four different angles, compared to the experimental data (Symbols 
as in Ref. The grey area is the region excluded from the fit (w, y/\t\ > 200 MeV). The dashed 
line is the 0{Q^) prediction. 

Thus the calculation of Ref. is not, strictly speaking, model independent: there is 
model input in the values used for Up and (3p. Using the same amplitude but fitting the 
polarizabilities, very good fits of the proton data in the low-energy regime (cu, \J\t\ < 200 
MeV) can be obtained. The central values for ap and (3p are similar to those employed in 
Ref. ^, but the uncertainty is larger: 



ar 



;i2.1±l.l)^°| X 10"^ fm^ 



(3p = (3.4± X 10'^ fm' 



(2) 



where statistical (l-c) errors are inside the brackets, and an estimate of the contribution 
from higher-order terms is given outside. These numbers are based on varying the upper 
bound on which data are fit from 160 MeV to 200 MeV, and on estimates of the 0{Q^) effect 
on 7p scattering. A sample of the best-fit results is shown, together with data, in Fig. ^ 
The fit has x^/d.o.f. = 170/131. These results are f ully compatible with other extractions, 
though the central value of f 5p i s somewhat higher 0[ . They are also consistent with 
the values predicted in Ref. |l6|, |l9|, where resonance saturation was used to estimate the 
0{Q'^) short-distance contributions. 
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The values (j21) can be compared with the recent re-evaluation of the Baldin Sum Rule by 
Olmos de Leon et al. which gives: 

ap + Pp = 13.8 ± 0.4. (3) 

This result overlaps the 1-a error ellipse of the fit ((21) , as shown in Fig. El Including the 
constraint Q in our fit leads to values for ap and f3p consistent with ((21) , but with smaller 
statistical errors, namely: 

ap = (11.0 ± 0.5 ± 0.2)+°| X 10~^fm^ 

(3p = (2.8±0.5t0.2)1[!:J X 10-Hm\ (4) 

In (0} we have left the systematic error unchanged, but have now included a second er- 
ror inside the brackets, whose source is the error on the sum rule evaluation ((SJ. The 
smaller overall statistical error in ((3} is achieved at the expense of additional (reasonable) 
assumptions about the high-energy behaviour of hadronic amplitudes, and also by using 
higher-energy data on the total proton photo-absorption cross section, since these are the 
two ingredients entering the Baldin Sum Rule result for ap + Pp. 

S. . , . . . , 

/3p 




FIG. 2: The 1-a error ellipse for our fit to the proton data with cj, < 200 MeV. Also shown 
is the band defined by the Baldin Sum Rule result ((21). The units of Qp and Pp are 10~^ fm^. 

The amplitude for Compton scattering on a nuclear target can also be calculated in the 
5FT, although the unraveling of scales is more subtle when more than one nucleon is present 
3|- In the single-nucleon sector a typical intermediate state has an energy denominator 
of 0{Q), but in multi-nucleon processes "reducible" intermediate states can have small 
energy denominators of 0{Q'^ /rriisf). The resulting infrared enhancement complicates the 
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perturbative expansion j2^. Furthermore, the very existence of nuclei imphes a breakdown 
of perturbation theory; hence the leading two-nucleon interactions should be summed to all 
orders, thereby building up a nuclear wave function In the^case of Compton scattering 
at energies of (9(m,r), this is the only resummation necessary 4], and the full amplitude, T, 
can be written in terms of a kernel K^^ which contains all irreducible (in the above sense) 
7A^A^ 'yNN graphs: 

T={i;\K.J^P). (5) 

The kernel can be calculated in xPT. Several other reactions involving deuterium have been 
successfully analyzed using analogous approaches ij]. However, for Compton scattering at 
energies of 0{m^/mi^) further resummations are necessary. In particular, such resumma- 



tions restore the Thomson limit for Compton scattering from deuterium 



In this 



energy regime it seems more appropriate to use a lower-energy EFT where pions have been 



integrated out j23, 



24| 



At momentum transfers of (9(m^), a formally-consistent power counting is emerging 
which organizes the nuclear interactions that give rise to the wave function This power 
count.ng . an in,p_t over Weinberg's ong.na, proposal Q, wMch ha. led to faWy 
accurate wave functions [1^. To the order we are working in the EFT, we require a deuteron 
wave function calculated to 0{Q'^) in the chiral expansion. We employ a wave function 
generated using the 0{Q^) xPT potential of Ref. j26|. (The A = 600 MeV wave function 
was chosen, but choosing A = 540 MeV instead produces very similar results.) In order 
to test the consistency of our error estimates, we have also computed results with other 
"realistic" wave functions, such as that obtained from the Nijm93 OBE potential ^22]. Below 
we will quote results only for the O^Q"^) wave function of Ref. j3] and for the Nijm93 wave 
function. The results for differential cross sections generated using obtained from other 
A^A^ potentials almost invariably lie between these two extremes. 

Meanwhile, the kernel K^^ is the sum of the single-nucleon Compton amplitude, in which 
one nucleon is a "spectator" to the Compton scattering, and "two-nucleon" contributions 
in which both nucleons are involved in the scattering of the photon. The amplitude for 
coherent Compton scattering on the deuteron was computed to 0{Q^) in Ref. j^. There are 
no free parameters to this order. The corresponding cross section is in good agreement with 
the Illinois data jol at 49 and 69 MeV, but under-predicts the SAL data 10 1 at 95 MeV. 



This calculation yields cross sections which agree well with the recent Lund data ^11,]. It 



also agrees qualitatively with potential-model calculations of 'jd scattering j28|, |29| . 

We have now extended our calculation to 0{Q^). The calculation contains both single- 
nucleon and two-nucleon contributions. For the first class of diagrams, we employ the 
single-nucleon amplitude of Ref. P]. The amplitude given there must be boosted from the 
Breit frame to the jd center-of-mass frame, which is straightforward once the jp amplitude 
has been decomposed into six invariant functions multiplying structures that transform in 
well-defined ways under boosts jj, |^ . 

Less straightforward are technical issues associated with the fashion in which nucleon 
recoil is included in the 'yd scattering calculation. These occur because the heavy-baryon 
formulation of xPT employed in this work expands about the limit of static nucleons, with 
nucleon recoil treated as a perturbative correction. The first issue has to do with the 
treatment of the very- low-energy region, uj ~ m^/M. At these energies the nucleons inside 
the deuteron are not static on the time-scales relevant to Compton scattering, and so a 
resummation of nucleon-recoil corrections must be performed. Below we will invoke this 
resummation when we calculate 'jd scattering for these low photon energies. The second 
issue has to do with ensuring that particle-production cuts appear at the correct position. 
In the absence of nucleon recoil the channel 'yd —>■ nd will open at a position which is in error 
by an amount / (2Mrf) . Since the opening of this channel can result in rapid variation of 
the Compton cross section with energy, this problem must be dealt with when we describe 
7P scattering for u close to m-Tr 0, II^- However, even at the highest energy for which 
Compton scattering on deuterium is calculated below, this small error in the position of the 
ird cut is a relatively minor effect. We emphasize that these subtleties are artifacts which 
arise due to the fact that nucleon recoil is treated as a perturbation in the 7N amplitude we 
have employed here P|. Such a treatment is in accord with the usual practice in the chiral 
EFT known as heavy-baryon xPT, but a perturbative treatment of nucleon recoil in not 
mandatory in the EFT. Thus these two problems in no way represent true limitations of a 
chiral EFT of Compton scattering. 

In addition, the two-nucleon diagrams at 0{Q'^) are added to the two-nucleon 0{Q^) 
diagrams shown in Ref. \^. These 0{Q^) graphs were generated by the leading chiral La- 
grangian. In going to 0{Q^) we include diagrams with one insertion from the sub-leading 
chiral Lagrangian. The coefficients of these vertices are essentially determined by relativis- 
tic invariance, and so these 0{Q^) two-body effects are suppressed by Q/rriN- No unknown 
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parameters occur in these graphs. However, two free parameters associated with the neu- 
tron polarizabihties do appear in the single-nucleon contribution: the shifts in the isoscalar 
polarizabihties from their 0{Q^) values. Indeed, of all the additional graphs which enter our 
new, 0{Q^), calculation only this effect from single-nucleon Compton scattering changes the 
cross section significantly. (Details can be found in Ref. 0.) 

We have fitted these two free parameters to the existing 'yd scattering data. As in the 
proton case, we impose a cutoff on the data that we fit when we extract a and (3. With 
a cutoff of uj, ^J\t\ < 200 MeV all of the 29 world data points 10, III are included. We 
also use a cutoff of 160 MeV, in which case all but the two backward-angle SAL points must 
be fitted. In both cases we float the experimental normalization for each experimental run 
within the quoted systematic error |8| , resulting in 22 (20) degrees of freedom for the fit P|. 

Using the wave function of Ref. Pq], we fit to data with u, < 160 MeV. This produces 
the isoscalar nucleon polarizabihties and the per degree of freedom given in the first line 
of Table U The is unacceptably large. This is driven mainly by the 49 MeV data from 
Illinois, and occurs because, as discussed above, for u ~ ml/rriN further resummations are 
necessary in the EFT. Here we adopt the strategy of including the dominant contributions 
that need to be added to our standard ^(Q^) result in this very-low-energy region in order to 
restore the Thomson limit for the 'yd amplitude [4] . This yields the results on the second line 
of Table IH An alternative strategy, namely dropping the 49 MeV data altogether, produces 
very similar central values for a a? and (3^- 
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(10-4 fm^) 


(10-4 fm^) 


xVd.o.f. 


Q' 


NLO xPT 


No 
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0.1 


2.36 




NLO xPT 


No 
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15.4 


-2.3 


1.95 


Q' 


NLO xPT 


Yes 


160 MeV 


13.0 


-1.8 


1.33 


Q' 


Nijm93 


Yes 


200 MeV 


16.9 


-2.7 


2.87 



TABLE L Results from different xPT extractions of isoscalar nucleon polarizabilities from yd 
scattering data. 



When the upper limit on which data is fitted is increased the central value oi — (5 



N 
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changes markedly, with a shghtly higher x^/d.o.f. (third hne of Table H}. Another key 
test involves examining the impact of the choice of deuteron wave function. The variability 
in the differential cross section due to the choice of wave function is ~ 10%. When the 
Nijm93 wave function is used the results are as shown in the fourth line of Table HI The 
high is again due to a failure to reproduce the 49 MeV data — this time even with the 
very-low-energy contributions included. While the Nijm93 wave function is not consistent 
with xPT, it does have the correct long-distance behavior. Such sensitivity to the choice of 
\'ip) is worrisome and merits further study. 

Putting these results together we conclude that our best fit, and error bars, for the 
isoscalar nucleon polarizabilities are: 

UN = (13.0± 1.9)1?| X 10~^fm^ 

= (-1.8 ± 1.9)12-1 X 10-^ fml (6) 

The statistical errors inside the brackets are obtained from the boundary of the 68% CL. 
region. The errors outside the brackets reflect the arbitrariness as to which data are included, 
and which deuteron wave function is employed. The results for the best-fit EFT are shown 
in Fig. El 

The results (EI) differ from those given in the published version of this work 18 ll . since 

6 p nnfT 

after publication we discovered two errors in our analysis of the 7d data |9|, |10|, llH- First, 
we mistakenly fitted the SAL data at a cm. energy of 95 MeV, whereas we should have 
employed a lab energy of 94.5 MeV jlQ]. (A similar error that ultimately has a much smaller 
impact was made with the data from Lund Q].) This causes a relatively small change in 
the results for most of our fits, but it does increase the central value for obtained with 
the Nijm93 deuteron wave function by 1.0 x 10~^ fm'^ over the value quoted in the published 



version of Ref. l31|| |32|. A larger effect occurs because there was an error in the program 

n 

used to generate the 0{Q^) results presented for 7d scattering in Ref. j3l|. This affects 
the 0{Q^) results by less than 1%. However, at 0((5^) the omission of this factor modifies 
the interference between the spin-dependent pieces of the single-nucleon amplitude and the 
0{Q^) two-body currents 2^. Correcting this mistake increases the predicted cross sections 



by about 10-15% with respect to those published in Ref. [3l|]. We emphasize that neither of 
these mistakes affects any of the results for 7p scattering reported here and in the published 
version of Ref. For more details see Ref. j^. 
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FIG. 3: Results of the 0{Q^) EFT best fit to the (lab and cm. as appropriate) differential cross 
sections for Compton scattering on deuterium at four different lab photon energies: 49, 55, 67, and 
94.5 MeV. The data are from Illinois 9] (circles), Lund 11] (diamonds) and SAL (squares) [lo| . 
The error bars represent the quoted statistical (only) uncertainties of these measurements. The 
solid line is the 0{Q^) calculation with on = 13.0 x 10~^ fm^, (3n = —1.8 x 10~^ fm^. The gray area 
is the region excluded from the fit (w, -^1*1 > 160 MeV). The dot-dashed line is the (parameter- free) 
0{Q^) calculation. 



Combining the revised numbers (jH)) with our results for and (3p we see that a wide range 
of neutron polarizabilities is consistent with a model-independent analysis of the current 
low-energy proton and coherent deuteron Compton data. Narrower ranges for the neutron 
polarizabilities can be obtained from this data at the expense of introducing model depen- 
dence. The values of a„ and /?„ extracted from data on the quasi-free process •yd — > •ynp 
using a theoretical model fall within this range However, there is clear statistical ev- 
idence in these data sets for q:„ being significantly larger than /3„ is consistent with 
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zero within our large error bars. Also, the results presented here provide no evidence for 
significant isovector components in a and (3. 

In conclusion, we have determined nucleon polarizabilities from a model-independent fit to 
low-energy Compton scattering on the proton and the deuteron. Our results are consistent, 
within error bars, with the recent extraction of oat ± from the Lund data using the 
detailed model of Levchuk and L'vov (But see also the values found using the data 

29|.). They are also consistent with the Baldin sum 
30|. The EFT can be improved by the introduction of 
an explicit A-isobar field, a complete treatment of the very-low-energy region, and a better 
understanding of the dependence of the results on the choice of deuteron wave function. An 
EFT study of the quasi-free deuteron process is also an important future step. 



of Refs. 0, Q and this model H 
rule results for + Pp and a„ + 
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